Models for Superfluid '^He in Aerogel 
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Several recent experiments find evidence of superfluidity of ^He in 98%-porous aerogel. The primary 
effect of the aerogel is that it scatters the quasiparticles of "^He. We find that many experimental 
findings are quantitatively understood by a relatively simple model that takes into account strong 
inhomogeneity of the scattering on a length scale of 100 nm. 
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The discovery of unconventional paring states in high- 
temperature superconductors has generated a lot of in- 
terest in impurity scattering in these materials. In par- 
ticular, the inhomogeneity of the scattering has been 
considered recently |l|]. However, both the experimen- 
tal and theoretical studies are difficult because of the 
complicated structure of these substances. Recently, a 
new possibility was opened for studying impurity effects 
on unconventional pairing states: superfluid "^He in very 
porous aerogel. This system has many advantages. For 
example, the pure state of superfluid '^He is absolutely 
pure in experiments, and it is theoretically very well un- 
derstood. A crucial parameter, the coherence length 
can easily be varied within a factor 5 by varying the pres- 
sure. The torsional oscillator experiments ||J^ and NMR 
experiments measure directly such basic quantities as 
the superfluid density, the pairing amplitude and the spin 
susceptibility. 

In this letter we give theoretical explanations for some 
of the experimental observations on superfluid ^He in 
aerogel. As a first attempt we study a model, where 
the aerogel is assumed as a homogeneous scatterer of the 
quasiparticles of '^He. This model gives predictions with 
a correct tendency, but it is insufhcient quantitatively. A 
"slab model" gives a clue that the inhomogeneity of the 
scattering is crucial for understanding the discrepancy. 
Based on that we construct a relatively simple model of 
inhomogeneous scattering that quantitatively can explain 
both the transition temperature and the pairing ampli- 
tude, and predicts an inhomogeneity length scale of 100 
nm. We also consider the upper limit for anisotropic 
scattering set by the NMR measurements. 

In the experiments the aerogel fills only 2% of the to- 
tal volume {V = 0.02), and its surface to volume ratio is 
A = 260,000 cm^^ Assuming naively that the ma- 
terial consists of a network of one-dimensional strands, 
we can from these numbers alone estimate the strand di- 
ameter AV/A — 3 nm. The distance between strands is 
V AttV I a = 20 nm. The mean free path for straight line 
trajectories is estimated as = 4/A = 150 nm. This is 
also the mean free path for quasiparticles of ■^He when 



the aerogel is filled with ■^He at millikelvin temperatures. 

Quasiclassical theory. — Because the volume fraction of 
the aerogel strands (including an inert layer of ^He atoms 
on the strands) is small, we neglect all effects that are 
linear in the volume fraction. In particular, we assume 
that the density, the Landau Fermi-liquid parameters, 
the coupling constant of the pairing interaction, and the 
dipole-dipole interaction constant are unchanged from 
the bulk. The changes of these parameters are of the 
same order of magnitude as the volume fraction because 
they arise from processes of relatively high energy and 
short length scale Much larger effects on superflu- 
idity arise from processes in the immediate vicinity of 
the Fermi surface. Scattering of quasiparticles from the 
aerogel strands modifies the superfluid state within the 
distance and causes an effect that is proportional to 
the ratio ^o/i, which approaches unity in 98%-porous 
aerogel. Here is the superfluid coherence length. It is 
deflned by = ^v-p /iiTk^TcQ, where Tco is the transition 
temperature in bulk ^He and up the Fermi velocity, 
is a function of pressure varying between 16 nm (melting 
pressure) and 77 nm (zero pressure). 

Because of the s-wave pairing of conventional super- 
conductors, their Tc and paring amplitude is nearly un- 
affected by non- magnetic scattering at ^o/^ ~ 1, only the 
Ginzburg-Landau coherence length gets shorter 0]. But 
in a p-wave superfluid like ^He, the scattering causes de- 
structive interference and leads to complete depression of 
superfluidity already at ^o/i ^ 1. 

All the models we discuss are quasiclassical. This 
means that the aerogel is modeled as a collection of in- 
coherent scattering centers at locations rj. Each center 
is assumed much smaller than but, similar to aerogel 
strands, they can be large in comparison to the Fermi 
wave length Ap = 2Tr/kp w 0.7 nm. For each scattering 
center, a fully quantum-mechanical treatment is allowed 
in principle, but we describe them phenomenologically 
by phase shifts and scattering cross sections. The in- 
terference of different scattering centers leads to weak 
localization corrections, which are small because aerogel 
has random structure and Xp/£ ^ 1. 
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FIG. 1. The transition temperature in aerogel relative to 
that in bulk, Tc/TcO- The horizontal axis is the coherence 
length ^0 = Tlvy I^Tik^Tco divided by L. The scale L is chosen 
so that the data sets coincide at the cross. The experimental 
results are from Refs. § (A, L = 36 nm), [|] (■, L = 25 nm), 
and (o, L = 24 nm). The lines correspond to the homoge- 
neous scattering model (HSM) to the slab model | fi5| , 
and to the isotropic inhomogeneous scattering model (IISM) 
with different sphere radii R and scattering profile parameters 
j- 



The coherence length is the only pressure dependent 
length scale in scattering models This implies that 
the calculated can be compared with experiments us- 
ing the scaling presented in Fig. 0. The vertical axis is 
the suppression of the transition temperature relative to 
the bulk, Tc/TcQ. The horizontal axis is divided by a 
length L. The scale L is a constant that characterizes 
each aerogel sample. By definition, L equals ^o(p) at the 
pressure p where Tc/Tco = 0.7. In other words, the hori- 
zontal scale is chosen so that all data sets coincide at the 
point (1.0,0.7). For the three different samples used in 
the experiments we find i = 36 nm L = 25 nm 
and L — 24: nm |^ . With this scaling the three data sets 
seem rather consistent with each other. 

In order to compare the amplitude A(r, r) of the order 
parameter, we study the suppression factor 0] 



(1) 



As before, the subscript refers to the bulk, i.e. to the 
case of pure '^He. The parameter t denotes the temper- 
ature relative to the transition temperature. An average 
over locations r is indicated by (...). 

We can construct a suppression factor Sp^ for the su- 
perfluid density ps in complete analogy with (^. How- 
ever, ps depends strongly on the Fermi-liquid parameter 
— 3(rncff / m — 1). Because the pressure dependence of 
spoils the scaling with discussed above, it is prefer- 
able to use the suppression factor Sp^ for the bare super- 
fiuid density p^ defined by ps = Ps/[l + ^FI{1 - Ps/p)], 
where p is the density of the liquid. 

The experimental suppression factors are plotted 
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FIG. 2. The suppression factors for the gap (5^2 ) and su- 
perfluid density {Sp^) as a function of squared Tc suppression, 
(Tc/Tcof. The upper two frames present the experimental 
data and the results of the homogeneous scattering model. 
The lower frames are results of the isotropic inhomogeneous 
scattering model. The numbers associated with curves and 
data points denote t, see equation (^. Born (sin^ 5o 0) 
and unitarity (sin^ So — 1) limits are shown by dashed lines 
in (a) and (b) at t — 1. All other curves are for the interme- 
diate case sin So = 0.5. 



against (Tc/Tco)^ in Fig. |(a, b). The NMR experiment 
measures Sa^ because the dipole-dipole interaction con- 
stant (7d 1^ is unchanged by scattering. Sp^ can be ex- 
tracted from torsional oscillator experiments. The t de- 
pendencies of S'a2 and Sp^ are qualitatively similar, but 
it is more pronounced in the latter. 

Homogeneous scattering model (HSM). — This is the 
simplest scattering model. The principal assumption is 
that the scattering probability is independent of the loca- 
tion. Additionally we assume that the scattering medium 
is isotropic, i.e. I is independent of the direction of quasi- 
particle momentum. These are just the standard assump- 
tions made in studying impurities in superconductors . 
We also neglect magnetic scattering because it does not 
seem to be important for the effects we consider. A con- 
venient property of the isotropic HSM is that both the 
Ginzburg-Landau (GL) theory and Leggett's theory of 
NMR [H have the same form as in pure ^He. Only the 
parameters of these theories have different values, as will 
be discussed below. 

The Ginzburg-Landau theory is formulated in terms of 
a free energy functional of the 3x3 matrix order param- 
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eter, A^i, where represents the spin components and i 
represents the orbital components of the pair state. The 
"bulk" terms are [|o|jl|] 

/bulk - aA*^^A^, + /3i|A^^A^,p + /32(A^«A;,)2 
'^P3A*^iA*^A^jAf^j + P4A*^^A^iAljAf^j 

+P,a;,a,,a,,a;^ . (2) 

The transition temperature Tc is determined by the con- 
dition a{Tc) = 0. Minimizing one finds the order 
parameter amplitudes A and free energies / of the vari- 
ous phases. For example, the polar, planar, and B phases 
have / = = -ka^/{4k(3i2 + 4/3345) with k = 1, 

2, and 3, respectively, where /3y... — /Si + l3j + and the 
A phase has / = aA^ — ~a^/ 4,^245- 
The coefficient a is given by |12] 



iV(0) 



-'^ cO -, 
n — 1 



, (3) 



where x = hvp/'^nTit,:, itr is the transport mean free 
path, and 2A^(0) the density of states at the Fermi sur- 
face. The suppression of Tc in the HSM is shown in Fig. 
|l| (^tr — S.7L). Its dependence on ^0 is the same as found 
for magnetic impurities in s-wave superconductors . 

For the coefficients Pi we make the additional assump- 
tion that only s-wave scattering is important, and obtain 
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Besides ^, h also depends on the s-wave scattering phase 
shift Sq. The effect of this fully quantum-mechanical de- 
gree of freedom on the suppression factors is shown by 
dashed lines in Fig. ||(a, b). However, calculations taking 
into account higher partial waves show that this depen- 
dence is essentially averaged out for large scatterers @ . 
The end result is approximately the same as if the phase 
shifts were random: sin^ (5o — > 0.5. Therefore we chose 
sin^ (5o = 0.5 in all other results of the HSM and IISM. 

The suppression factors of the HSM are essentially the 
same for different superfluid phases. The difference be- 
tween A and B phases is negligible also in Sp^ , when the 
average of the tensor *ps is used for the A phase. The 
stability of A and B phases depends on strong coupling 
corrections A/?^'^, which are not known. Assuming they 
remain constants, the B phase is favored by increasing 
scattering. No new phases are stabilized. 



We conclude the HSM by noting that it works in the 
right direction for all Tc, A, and ps, but quantitatively, on 
the level of accuracy we are accustomed to in superfluid 
■^He, it is clearly inadequate. 

Slab model. — This model considers ^He in a gap of 
thickness D between two diffusely scattering planes. The 
dashed line in Fig. |^ shows Tc calculated in Ref. p^ ] 
(D = 2.95T). The agreement with measurement is much 
better than for the HSM. In particular, Tc suppression 
is quadratic at small ^0 compared to linear in the HSM. 
Generally, this feature arises from regions that have no 
scattering nearby, such as the center of the slab. The 
suppression of (A^) also is in better agreement with ex- 
periments than in the HSM . The principal deficiency 
of the slab model is its strong anisotropy. We estimate 
(see below) that in order to be in agreement with the 
measured NMR shift, the normal direction of the slab 
has to vary randomly on a length scale that is smaller 
than the thickness D. This contradiction prompts us to 
look for a better model. 

Isotropic inhomogeneous scattering model (IISM). — 
Experimentally the suppression of Tc seems quadratic 
at small ^0 (Fig- 0)- This implies that real aerogel has 
voids, i.e. regions of negligible scattering. In order to 
construct a model that is feasible in calculations, we 
make two basic simplifications, (i) Instead of a random 
distribution of voids, we consider a periodic lattice of 
them, (ii) The unit cell of this lattice is approximated 
by a sphere. In more detail, the boundary condition 
is that a quasiparticle escaping from the sphere will be 
returned there at the diametrically opposite point but 
its momentum is unchanged. (A phase shift similar as 
in Bloch wave functions should be added in the case of 
nonconstant phase.) In addition to the radius R of the 
sphere, we need to specify how the density n{r) of scat- 
tering centers is distributed in the sphere. When n(r) 
depends only on the radial coordinate, i.e. n(r) = n(r), 
the model is completely isotropic. We study polynomial 
forms n{r) = c\{r / RV — i{r / Ry^"^ / {j + 2)]. We show re- 
sults in Figs. |l] and ||(c, d) for a steep (j = 8) and a slow 
(j = 2) profile. The suppression factors are for an in- 
homogeneously distorted B phase, but extrapolating the 
experience from two previous models, the A-type phase 
would be quite similar. 

In spite of the inhomogeneity, there is a single well de- 
fined Tc at which the order parameter becomes nonzero. 
Due to the proximity effect, Tc is determined collec- 
tively by the whole sample, although the weight of high- 
scattering regions far from low-scattering ones is expo- 
nentially small. Anyway, the transition can be described 
as "broadened" if (A'^)(T) is strongly nonlinear below 
Tc. This is the case for a slow profile [dashed lines in 
Fig. |(d)] and for large Tc/Tcq: in Fig. |(c, d) this shows 
up as strong t dependence of S\2 . In contrast, ( A^) (T) is 
nearly linear over the whole temperature range for small 
Tc/TcQ. [In this case, the t dependence of arises 
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mostly from nonlinearity of the reference Aq(T).] 

The IISM reduces to the HSM in the hmit of small R. 
This means that the true distribution of the scattering 
centers is irrelevant as long as the average scattering over 
a length scale remains the same 

We see that the IISM is in much better agreement with 
experiments than the HSM when 5L and j « 8. The 
magnitude and t dependence of Sa^ and most of Tc(^o) is 
well accounted for. There is a small systematic deviation 
that experimentally both and drop more rapidly 
with increasing than in the model. We believe that this 
difference arises from the periodicity assumption in the 
IISM. In real aerogel there are fluctuations of all length 
scales, and with increasing the length scale of most 
relevant fluctuations also increases. This is consistent 
with the observed deviations which imply an increasing 
effective R for increasing • 

The Ps{T) measurement [Fig. ^(b)] shows consider- 
ably more nonlinearity than the NMR measurement [Fig. 
11(a)] . A possible explanation for this is that the struc- 
ture of the aerogel is different in the two experiments, the 
former corresponding to a smaller j. In order to confirm 
such a hypothesis, both samples should be studied with 
the same measuring technique. 

The large scattering fluctuations predicted by the IISM 
seem to us natural and essentially unique explanation 
of the measured suppression factors. However, the pre- 
dicted effective void radius « 0.8R ~ 100. ..150 nm is very 
large compared to the estimates based on small angle x- 
ray scattering |^. This problem remains open. 

Anisotropic HSM. — According to our introductory es- 
timate, the aerogel consists of randomly oriented strands 
of diameter 3 nm and length La ~ 20 nm. This kind of 
anisotropy can have several consequences on the super- 
fluid state. It is known that strongly anisotropic scat- 
tering, such as in the slab model, can stabilize the A 
phase, and one can ask if the aerogel strands could do 
the same. We have studied this in the limit La ^ foi 
where we recover the GL functional (||) of the isotropic 
HSM since the anisotropy is averaged out on the scale fo- 
However, the coefficients /3i are modified. We find that 
anisotropic backscattering, preferentially perpendicular 
to the strands, can stabilize the A-phase low pressures, 
where the B-phase is otherwise stable p^ . 

The anisotropy couples to the orbital part of the order 
parameter. For example, the a-term in the GL functional 
is modified to aifcA*^^^^. Let us study this in the A 
phase, where the orbital and spin parts are described by 
1 and d vectors, respectively. There are two possibilities 
[ p^ . (i) For weak anisotropy the dipole-dipole coupling 
between 1 and d keeps them aligned to each other, (ii) 
For strong anisotropy, 1 is driven to vary randomly on a 
scale Lo where d is still nearly constant. Using similar 
estimates than Imry and Ma [|l^ we find that NMR fre- 
quency shift in the former state is unchanged relative to 



HSM, but it is reduced to essentially zero in the latter. 
Experiments clearly point to the former state 0|. Our 
scattering estimate also favors this state, but the margin 
is rather small: if La « 50 instead of 20 nm, the latter 
state would be favored. The proximity of the transition 
gives natural explanation to the observed sudden extinc- 
tion of the NMR shift as a function tipping angle [^,|l^ . 

Conclusions. — Superfiuid "^He in aerogel is in many re- 
spects an ideal system to study impurity effects in un- 
conventional superfluidity. We find, in particular, that 
the standard impurity model is robust in the sense that 
large fluctuations in the scattering are needed in order to 
get such substantial deviations as seen experimentally. 

We thank R. Hanninen and T. Setala for help in nu- 
merical calculations. 
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